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RA Z MVE MSHI MA SIT A Sie 


郭 江 鹏 , EMS, SOR, H Mk, 李建平 , WAR, BAM 
(河北 农业 大 学 机 电工 程 学 院 , 河 北 保定 071000) 


摘 要 : 针对 果园 多 风 道 喷雾 机 内 部 气流 分 布 不 均 导 致 由 出 风口 吹出 的 气流 北 乱 、 影 响 使 雾 滴 在 果树 冠 层 
上 均匀 沉积 的 问题 ， 对 多 风 道 喷雾 机 内 部 导 流 板 长 度 参数 进行 了 优化 。 应 用 计算 流体 动力 学 (Computational 
Fluid Dynamics, CFD) 技术 ,基于 Star-CCM+ 软 件 对 喷雾 机 送 风 系统 内 部 气流 进行 了 模拟 分 析 ， 得 到 出 风口 
1~6 的 风速 在 不 同 导 流 板 长 度 的 标准 差分 别 为 0.7468、0.6776、1.4441、5.1305、4.5768 和 0.8209。 对 风速 标 
准 差 较 大 的 出 风口 3、 出 风口 4、 出 风口 5 进行 响应 面 分 析 ， 最 终 确定 导 流 板 1 长 度 200.00 mm、 导 流 板 2 长 
度 60.00 mm、 导 流 板 3 长 度 50.00 mm 为 最 优 参 数组 合 。 在 最 优 组 合 参数 下 ， 计 算得 到 对 称 出 风口 3 和 出 风口 
6 的 风速 值 分 别 为 39.135 和 41.320 m/s， 相 对 偏差 为 5.58%; 出 风口 4 和 出 风口 5 的 风速 值 分 别 为 33.022 和 
34.328 ms， 相 对 偏差 为 3.95% ， 符 合 设计 要 求 。 室 内 风速 试验 结果 表明 ， 在 距离 喷雾 机 出 风口 1.25 m kh, 
风 场 风速 由 上 层 到 下 层 逐 渐 增 大 ， 实 现 风 场 按 果树 冠 层 形状 分 布 ， 喷雾 机 左右 两 侧 风 场 对 称 分 布 ， 气 流 分 
布 均 匀 。 果 园 多 风 道 喷雾 机 设计 满足 要 求 ， 可 为 同类 设计 提供 参考 。 
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1 引 言 将 模拟 与 试验 相 结合 进行 探索 是 应 对 这 些 影响 因 

素 的 有 效 方法 。 在 研究 过 程 中 ， 计 算 流体 力学 

l a 

素 的 影响 ， 如 农药 剂量 和 喷雾 量 、 药 液 分 布 、 风 ee 
量 、 气 流 模式 、 喷 雾 机 速度 、 气 象 条 件 以 及 农 作 “了 充分 的 应 用 ， 用 以 帮助 解决 复杂 的 过 程 。 

物 特性 等 5 ， 需 最 大 程度 地 提高 生物 学 功效 、 HAY A R E AH CFD 技术 广泛 应 用 在 风 送 


果园 喷雾 作业 过 程 复杂 ， 喷 雾 效 果 受 很 多 因 


最 小 化 残留 物 ， 并 与 环境 污染 保持 平衡 54。 果 喷雾 机 的 设计 研究 上 。Wang 等 ”通过 对 轴 流 风 
园 喷 雾 涉 及 大 量 的 喷雾 器 参数 设置 和 环境 参数 ， ”机 模拟 数值 分 析 得 出 当 进 口 速度 在 25 一 55 m/s, 
仅 使 用 试验 方法 很 难 进行 分 析 且 成 本 昂贵 9" ， 扩散 角 在 4.5? 或 5. 时， 压力 损失 最 小 ， 而 随 着 
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A ARE SS, KIRKEE. EB 
AEU kih T — ak A E RR pi Ze rl AIK BSE BL 
通过 理论 设计 和 虚拟 样机 技术 完成 整 机 结构 设计 
和 各 关键 部 件 参 数 的 确定 ， 并 通过 田间 试验 验 
证 性 能 。 周 良 富 等 ' 设计 了 组 合 圆 盘 式 果园 风 
送 喷雾 机 ， 并 设计 田间 试验 得 到 了 最 佳作 业 参 
Bl. RAB! 利用 CFD 技术 ， 针 对 离心 风机 出 风 
腔 进行 风 场 分 析 ， 结 果 显 示 在 风机 腔 中 加 装 带 圆 
角 的 导 流 块 ， 可 以 使 风机 风量 提高 0.27% 一 
3.69%。 丁 天 航 等 "针对 单 风 机 流 道 果 园 喷雾 机 
气流 场 分 布 不 均 的 问题 ， 采 用 CFD 数值 模拟 ， 
对 比分 析 了 单 双 风 道 的 气流 场 。Endalew 等 gÈ 
立 了 一 种 基于 CFD 的 果园 风 送 模型 ， 该 模型 能 
人 够 直接 捕获 果树 冠 层 形状 及 其 对 自然 风 和 喷雾 气 
流 的 局 部 影响 ， 从 而 得 出 更 为 真实 的 结果 。 
Nuyttens 等 “通过 CFD 技术 建立 喷雾 机 三 维 雾 
WRR, ， 结 合 果树 冠 层 及 环境 因素 对 喷雾 机 
风 场 进行 模拟 分 析 。Badules 等 "建立 了 一 种 
CEFD 模 型 ， 用 于 研究 喷雾 需 作 业 速 度 对 横流 式 风 
送 喷雾 机 外 部 三 维 气流 的 影响 。Hong 等 ”开发 
了 一 个 集成 的 CFD 模型 用 于 预测 由 风 送 喷雾 机 
吹出 的 气流 在 树木 冠 层 内 部 和 周围 的 速度 分 布 ， 
该 模型 通过 滑动 网 格 技术 模拟 喷雾 需 的 运动 ， 并 
且 在 计算 域 中 将 树木 冠 层 定义 为 虚拟 多 孔 介 质 ， 
利用 该 模型 分 析 有 树枝 的 树木 冠 层 内 部 气流 速度 
分 布 。Zheng 等 "利用 CFD 技术 对 喷雾 机 外 部 
风 场 进行 模拟 分 析 ， 在 实验 室 条 件 下 进行 验证 并 
测试 。Holownicki 等 "开发 了 一 个 集成 的 CFD 


的 分 布 ， 而 多 风 道 风 送 式 喷雾 机 的 内 部 导 流 结构 
布置 不 合理 会 直接 导致 出 风口 风速 分 布 不 均匀 的 
现象 。 为 对 果园 多 风 道 喷雾 机 内 部 风 道 进行 优化 
设计 ， 得 到 内 部 导 流 板 最 佳 参数 组 合 ， 本 研究 应 
用 CFD 技 术 ， 基 于 StarCCM+ 软 件 对 果园 多 风 道 
喷 筋 机 送 风 系 统 内 部 气流 进行 模拟 分 析 ， 采 用 数 
值 模拟 分 析 的 方法 研究 送 风 系 统 内 部 气流 场 分 布 
规律 ， 对 各 个 参数 下 的 简化 模型 进行 数值 模拟 分 
析 ， 后 期 观察 送 风 系 统 内 部 流 场 和 出 风口 风速 ， 
最 终 确 定 风 道内 部 导 流 板 长 度 ， 实 现 喷 筋 机 两 侧 
风 道 出 风 对 称 分 布 。 


2.1 风 道 设计 


多 风 道 出 口 装置 是 连接 涡轮 式 离心 风机 的 关 
键 零 部 件 ， 将 离心 风机 产生 的 气流 分 配 到 各 风 道 
出 口 。 传 统 的 风 送 式 果 园 喷 雾 机 多 由 轴 流 风机 或 
离心 风机 产生 喷雾 所 需 气流 ， 通 过 导 风 装置 〈 如 
导 流 板 ) 改变 气流 场 的 分 布 。 此 类 导 风 装置 的 出 
口 面积 较 大 ， 从 而 导致 气流 场 的 分 布 主要 呈 半 圆 
形 扩散 ， 部 分 由 风机 产生 的 气流 不 能 有 效 到 达 葛 
标 区 域 ， 造 成 气流 损失 ， 进 而 导致 出 风口 处 的 气 
流风 速 降低 。 本 研究 设计 了 多 风 道 出 口 装置 ， 由 
多 个 相同 截面 的 内 部 风 道 组 成 ， 保 证 气流 能 够 精 
准 有 效 地 流向 靶 标 所 在 的 方向 。 根 据 前 期 对 果树 
冠 层 的 测量 数据 ， 对 风 道 分 布 进行 计算 。 将 果树 
和 冠 层 近似 看 作 圆 锥 形状 ， 对 冠 层 进行 分 层 处 理 ， 
上 中 下 层 冠 层 半 径 分 别 为 0.40、0.80 和 1.10 m, 


模型 用 来 预测 风 送 式 喷 雾 器 吹出 的 树冠 内 和 周转 
气流 的 速度 分 布 ， 采 用 滑 移 网 格 技术 模拟 喷头 运 
动 ， 将 树冠 模型 定义 为 虚拟 多 孔 介质 ， 通 过 将 
CFD 模 拟 结 果 与 之 前 的 测量 结果 进行 比较 ， 验 证 
CFD 模 型 。 翟 长 远 等 "应 用 CFD 流体 仿真 技 
术 ， 对 塔 式 喷雾 机 出 风口 风 场 进行 建 模 与 试验 ， 
研究 了 汕 流 模型 和 计算 域 尺寸 对 风 场 分 布 结果 的 
影响 特性 。 

以 上 研究 针对 不 同类 型 喷雾 机 进行 的 CFD 
风 场 研究 发 现 ， 模 型 构建 的 不 同 直接 影响 着 风 场 


依据 公式 (1) 对 上 中 下 冠 层 体积 进行 计算 。 

V = smh (R tre+ Rr) (1) 

其 中 ， 大 表示 冠 层 体积 m; 7 表示 上 中 
下 层 之 间 的 高 度 差 , m， 取 7=1.00m ( 树 高 在 
3.50 m 左 右 ， 距 地 面 0.50 m 处 为 果树 侧枝 生长 初 
始点 ， 因 此 剩余 3.00m 上、 中 、 下 层 之 间 的 平均 
高 度 h 取 1.0m); RÆ; TRÍ, m; 7 表示 上 
底 半 径 ，m。 

经 计算 得 到 果树 下 、 中 、 上 层 的 冠 层 体积 分 
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别 为 2.8、1.1 和 0.3m， 在 设计 出 风口 时 ， 下 层 
两 个 出 风口 需 占 到 整个 圆 弧 面 的 1/12， 中 层 和 上 
层 出 风 口 由 圆 弧 面向 上 延伸 ， 进 行 4 等 均 分 ， 由 
内 部 三 个 导 流 板 控 制 风 量 ,， 根 据 第 3 节 送 风 系 统 
内 部 气流 场 模拟 分 析 确 定 由 左 到 右 的 三 块 导 流 板 
的 长 度 分 别 为 200.00、60.00 和 50.00 mm， 如 图 1 
所 示 。 设 计 内 部 风 道 时 ， 为 避免 气流 在 较 大 转弯 
处 产生 汕 流 现象 ， 设 置 通用 导 流 片 ， 保 证 气流 能 
够 平稳 通过 ， 进 而 使 气流 速度 损失 最 小 。 

查阅 相关 资料 ， 结 合 实 地 调研 河北 省 保定 市 
曲阳 县 河北 绿 阳 现代 农业 园区 (38°N, 114° E) 
内 的 苹果 园 ， 芋 果树 行距 为 3.30 m。 喷 雾 机 工作 
最 佳 喷 雾 距离 为 1.25 m 2 ， 再 通过 农机 农艺 相 
结合 ， 设 计 得 到 多 风 道 喷雾 机 风 道 部 分 总 体 尺 十 
宽 1.00 m， 高 2.00 m， 利 用 AIP (Autodesk In- 
ventor Professional) 三 维 软件 进行 建 模 ， 得 到 风 
道 主体 部 分 结构 示意 图 ， 如 图 1 所 示 。 


Dn vA Bw 


注 ; 1. 出 口 1; 2. 出 口 2; 3. 出 口 3; 4. 出 口 4; 5. 出 口 5; 6. 出 口 6; 
7. 导 流 板 1; 8. 导 流 板 2; 9. 导 流 板 3 

Al 多 风 道 喷雾 机 六 风 道 主体 结构 示意 图 
Fig. 1 Schematic diagram of main structure of six air ducts 


multi-air-sprayer 
2.2 出 风口 设计 


果园 多 风 道 喷雾 机 利用 风机 产生 的 高 速 气流 
对 雾 滴 进 行 二 次 雾 化 ， 使 得 雾 滴 能 够 与 靶 标 充分 
接触 ， 故 风机 产生 的 气流 应 与 雾 滴 充分 接触 。 由 
于 圆 形 风 道 出 风 为 直线 ， 无 法 产生 射流 扩散 ， 根 
据 图 2 所 示 的 喷雾 示意 图 ， 设 计 了 一 种 鸭 嘴 型 组 
合 喷头 出 风口 ， 出 风口 长 度 为 360 mm， 宽 度 为 
45 mmm， 扇面 扩散 角 为 55.0*， 以 实现 气流 在 果树 


冠 层 上 的 全 和 覆盖。 该 设计 包括 进 风 口 、 鸭 嘴 型 出 
风口 、 喷 管 、 喷 头 、 固 定 片 等 组 成 ， 结 构 框 图 如 
图 3 所 示 。 


注 : 1. 鸭 嘴 形 喷雾 装置 ; 2. 延 伸 简 ; 3. 六 出 口 风 道 ; 4. 风 机 


图 2 多 风 道 喷雾 机 喷雾 示意 图 


Fig. 2 Schematic diagram of multi-air-sprayer Spray 


注 : 1. 喷 管 ; 2. 喷 头 ; 3. 鸭 嘴 型 出 风口 ; 4. 固 定 片 ; 5. 进 风口 
图 3 鸭 嘴 型 组 合 喷头 出 风口 


Fig. 3 Outlet of duckbill type combined nozzle 

在 本 设计 中 ,由 和 矩形 截面 的 内 部 风 道 转 为 圆 
形 截面 的 送 风 通 道 ， 再 由 圆 形 截面 到 鸭 嘴 型 出 风 
口 的 过 程 中 ， 可 以 有 效 增 大 风 压 并 扩大 出 风口 的 
喷 筋 幅 宽 ， 使 得 喷头 筋 滴 与 蝎 嘴 型 出 风口 的 气流 
能 够 充分 接触 且 气 流 分 布 均匀 ， 二 次 雾 化 效 
果 明 显 ， 气 流 携带 雾 滴 均 匀 喷 洒 到 靶 标 区 域 ， 同 
时 翻转 叶片 ， 使 果树 叶片 正 、 反 面 均匀 着 药 。 


3 内 部 气流 场 模拟 分 析 


3.1 影响 因素 分 析 


果园 多 风 道 风 送 式 喷 雾 机 中 上 层 出 风 口内 部 
风 道 由 三 块 不 同 长 度 的 导 流 板 分 割 开 ， 三 块 导 流 
板 能 够 调节 离心 风机 产生 的 高 速 气流 ， 实 现 对 喷 
雾 机 内 部 送 风 系 统 气 流 场 的 分 配 ， 如 图 1 所 示 。 
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离心 风机 产生 的 高 速 气流 在 喷雾 机 风 道 的 导 流 作 
用 下 被 定向 输送 至 喷雾 机 各 出 风口 ， 导 流 板 的 导 
向 过 程 在 很 大 程度 上 影响 着 喷雾 机 的 作业 效果 。 
不 同 的 导 流 板 长 度 会 影响 着 喷雾 机 内 部 气流 的 分 
布 ， 产 生 不 同 的 喷 筋 效果 ， 同 时 影响 着 喷雾 机 的 
作业 效率 。 

经 分 析 ， 影 响 喷 筋 机 内 部 气流 的 主要 因素 为 
三 块 导 流 板 的 长 度 ， 其 中 ， 导 流 板 1 的 长 度 直接 
影响 着 出 口 3、 出 口 4 以 及 出 口 5 的 风速 及 风量 ， 
导 流 板 2 和 导 流 板 3 则 对 出 口 3、 出 口 4 及 出 口 S 
处 的 气流 进行 导向 ， 如 图 1 所 示 。 为 确定 最 佳 组 
a, 设置 3 因素 3 水 平 试 验 水 平 表 ， 导 流 板 长 度 
需 满 足 系统 对 风量 的 要 求 ， 试 验 因素 下 限 由 风量 
要 求 确定 ， 该 系统 风量 大 小 应 满足 3.61 mYs。 按 
照 不 同 导 流 板 规律 依次 递减 ， 设 计 三 块 导 流 板 
长 度 初 始 长 度 分 别 为 : 导 流 板 1 的 长 度 工 为 
225.00 mm， 导 流 板 2 的 长 度 疡 为 85.00 mm， 导 
TRA 3 WK EEL, 65.00 mm， 如 表 1 所 示 。 使 用 
Design-expert 软件 设计 回归 试验 ， 试 验方 案 表 如 
表 2 所 示 。 
表 1 多 风 道 喷雾 机 送 风 系统 模拟 试验 的 因素 水 平 表 
Table 1 Factor table of simulation test for air supply system 


of multi-air-duct sprayer 


me 试验 因素 
L /mm L,/mm L,/mm 
1 225.00 85.00 65.00 
0 200.00 60.00 50.00 
-1 175.00 35.00 35.00 

表 2 试验 方案 设计 表 

Table 2 Test scheme design 

试验 序号 L /mm L,/mm L /mm 
1 175.00 35.00 35.00 
2 200.00 60.00 35.00 
3 225.00 35.00 50.00 
4 175.00 60.00 35.00 
5 200.00 60.00 50.00 
6 225.00 60.00 35.00 
这 175.00 60.00 65.00 
8 200.00 85.00 35.00 
9 225.00 60.00 65.00 


3.2 模型 建立 与 仿真 试验 


利用 AP 软件 对 喷雾 机 模型 进行 初始 建 模 ， 
并 按 试验 要 求 更 改 喷雾 机 导 流 板 长 度 ， 将 建 好 的 
13 组 流体 域 模型 依次 导入 Star-CCM+ 软 件 中 。 在 
os # 行 表面 修复 、 布 尔 运 
算 、 边 界 条 件 设 定 、 物 理 模型 定义 、 网 格 划 分 以 
eam ed 
喷雾 机 流体 域 模型 (图 4) 采用 多 面体 网 格 ， 并 
对 交界 面 和 风机 壁面 进行 局 部 加 密 处 理 。 在 模型 
设置 完成 后 ， ee 
(图 4) ， 并 根据 各 监测 点 建立 速度 监测 报告 ， 和 
出 各 出 风口 处 的 平均 气流 风速 计算 结果 ， 同 时 选 
取 中 间 截 面 用 以 观察 后 处 理 中 速度 矢量 分 布 图 。 


oint6 了 Points Point4 Point3 


注 : Point 1 为 监测 点 1, 以 此 类 推 
图 4 多 风 道 喷雾 机 流体 域 模 型 各 监测 点 分 布 图 


Fig. 4 Distribution charts of monitoring points for multi-air- 


duct sprayer fluid domain model 


4 仿真 结果 与 分 析 


4.1 基于 Star-CCM+ 的 内 部 流 场 仿真 分 析 


模拟 仿真 结 ， 将 Star-CCM+ 软 件 中 得 到 
的 各 监测 点 的 "o 报告 以 及 各 模型 截面 的 速 
度 矢量 图 导出 ， 如 图 5。 将 仿真 计算 的 结果 进行 
汇总 并 分 析 ， 图 5 中 能 明显 看 出 风 道内 部 气流 变 


化 情况 。 
整理 汇总 分 析 各 监测 点 的 速度 监测 报告 ， 结 
果 如 图 6。6 组 监测 点 按照 1~6 顺 序 ， 在 不 同 导 


流 板 长 度 的 速度 标准 差分 别 为 0.7468、0.6776、 
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EE ms) EE =) 速度 (ms 
191.749 pas 194.925 
73.399 mas 75.940 
55.049 56.606 56.955 
36.700 pRB? 37.970 
18.350 i 18.985 
0.000 pooo 10.000 
(a)175+35+50 (b)200+35+35 (¢)225+35+50 
速度 /(m- s7) 速度 ((m- s7) 速度 /(m- s 
85.979 94.270 94.667 
38.783 75.416 75.733 
51.588 56.562 56.800 
34.392 37.708 37.867 
17.196 18.854 18.933 
10.000 10.000 10.000 
(d)175+60+35 (e)200+60+50 (j)225+60+35 
AE s7) ER (ms) 速度 (ms 
pane 94.821 
73.696 75.152 75.857 
55.272 56.364 56.893 
36.848 37.576 37.928 
18.424 18.788 18.964 
0.000 10.000 10.000 
(h)175+60+65 (4)200+85+35 (j)225+60+65 


注 : 以 导 流 板 1+ 导 流 板 2+ 导 流 板 3 作为 参数 组 合 编码 规则 ,例如 200+65+35 表 示 导 流 板 1 长 度 为 200.00 mm, Fi 2 KZA 65.00 mm, 


导 流 板 3 长 度 为 35.00 mm 


AS 多 风 道 喷雾 机 模型 截面 速度 矢量 分 布 图 


Fig. 5 Velocity vector distribution of model section multi-air-duct sprayer 


1.4441, 5.1305, 4.5768 和 0.8209。 其 中 Pointi, 
Point2 、Point6 这 三 个 监测 点 的 风速 值 标准 差 均 
小 于 1， 表 明 导 流 板 长 度 的 变化 对 出 风口 1、 出 
风口 2、 出 风口 6 的 速度 变化 影响 较 小 可 忽略 不 
计 ， 而 Point3 、Point4 、Points 这 三 个 监测 点 的 风 
速 值 标准 差 较 大 ， 表 明 随 着 导 流 板 长 度 的 变化 ， 
出 风口 3、 出 风口 4、 出 风口 5 处 的 风速 受 影响 程 
度 大 。 


结果 分 析 与 验证 


为 得 到 导 流 板 长 度 最 优 参数 组 合 ， 采 用 响应 
面 法 (Response Surface Methodology, RSM) 对 
Point3 、Point4 Points 这 三 个 点 的 平均 模拟 风速 
值 进行 进一步 研究 ， 将 三 个 点 的 风速 响应 值 记 为 
了 五、 了 对、 了 影响 因素 导 流 板 1 的 长 度 记 为 A、 导 


a 


A 
Pointl Point2 Point3 Point4 i Point6 


出 风口 监测 点 
图 6 喷雾 机 内 部 风 道 流体 域 模型 风速 计算 结果 


Fig. 6 Fluid velocity calculation results in internal duct of 


the sprayer 


流 板 2 的 长 度 记 为 B、 导 流 板 3 的 长 度 记 为 C， 进 


行 回 归 分 析 ， 得 出 回归 方程 。 
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42.1 出 风口 3 
对 出 风口 3 处 的 平均 模拟 风速 值 进行 分 析 ， 
得 到 回归 方程 如 公式 (2)。 
Y,=39.13-1.07 X A-0.13 X B-0.17 X C+ 
0.13 X AB +0.19 X AC+0.16 X BC- 
1.84 X A?- 1.34 X B?-0.58 XC (2) 


回归 方程 的 显著 性 及 方差 分 析 如 表 3 所 示 。 
模型 P 值 ~<0.01， 为 极 显 著 ， 决定 系数 R= 
0.9979， 表 明 模 型 拟 合 效果 较 好 ， 可 用 来 分 析 及 


表 3 多 风 道 喷雾 机 出 风口 3 平均 风速 值 响应 面 回归 方程 
的 显著 性 及 方差 分 析 结 果 
Table 3 Significance and variance analysis of response sur- 


face regression equation for outlet 3 of multi-air-duct sprayer 


方差 来 源 ”平方 和 自由 度 HF FE  P 值 显著 性 


模型 35.38 3.93 362.49 <0.0001 显著 
4- 导 流 板 1 长 度 9.24 9.24 851.81 <0.0001 极 显著 
B- 导 流 板 2 长度 0.13 0.13 12.36 0.0098 极 显著 
C- 导 流 板 3 长 度 0.24 0.24 22.53 0.0021 极 显著 


AB 0.067 0.07 6.11 0.0427 显著 
AC 0.14 0.14 13.26 0.0083 极 显著 
BC 0.11 0.11 9.72 0.0169 ”显著 


B 7.53 7.53 694.51 <0.0001 极 显著 
C 1.42 1.42 131.36 <0.0001 极 显著 
残 差 0.08 0.01 
失 拟 项 0.08 0.03 
纯 误 差 0.00 
总 离 差 35.45 16 


9 
1 
1 
1 
1 
1 
1 
a 14.32 1 14.32 1320.81 <0.0001 极 显 著 
1 
1 
7 
3 
4 


0.00 


ms 


模拟 风速 值 / (mr 
模拟 风速 值 / G 


w 
A: FALHI /mm 


3s yé 


(a) 导 流 板 1 与 导 流 板 2 交互 响应 面 


(b) 导 流 板 1 与 导 流 板 3 交 互 响 应 面 


预测 响应 值 。 该 模型 的 一 次 项 4、B、C 极 显著 ; 
3 个 因素 之 间 两 两 交互 4C 极 显著 ，4B、BC 显 
著 ; 二 次 项 4、PP、C 极 显著 。 出 风口 3 处 的 平 
均 风 速 值 头 作 为 评价 标准 ， 在 数值 分 析 中 作为 优 
(Ate, WEYER, ARES RI, Bi 
板 2、 导 流 板 3 的 长 度 ， 数 学 模型 表示 如 
公式 (3)。 

max: Y, =f (A,B,C) 

35<C< 60 

为 检验 试验 结果 的 可 靠 性 ， 采 用 最 优 导 流 
板 参数 组 合 进行 3 次 平行 仿真 试验 ， 得 到 的 出 
风口 3 平均 风速 值 分 别 为 38.789、39.420 和 
38.517 m/s， 平 均值 为 38.909 ms， 结 果 与 预测 值 
接近 ， 相 对 误差 为 0.58%， 表 明 该 模型 能 够 准确 
的 预测 仿真 结果 。 

出 风口 3 平均 风速 值 的 响应 面 结 果 如 图 7 所 
示 。 由 图 7 可 以 直观 地 发 现 导 流 板 1、 导 流 板 2、 
导 流 板 3 之 间 的 交互 作用 ， 两 个 因素 同时 作用 下 
风速 变化 值 明显 ， 交 互 作用 明显 。 分 析出 风口 3 
处 的 平均 风速 值 ， 在 各 试验 因素 的 取 值 范围 内 ， 
求 得 最 优 的 导 流 板 长 度 组 合 为 导 流 板 1 长 度 
200.00 mm、 导 流 板 2 长 度 60.00 mm、 导 流 板 3 
K Æ 50.00 mm， 出 风口 3 平均 风速 值 为 
39.135 m/s. 


(Cc) 导 流 板 2 与 导 流 板 3 交互 响应 面 


图 7 多 风 道 喷雾 机 出 风口 3 平均 风速 值 的 响应 面 


Fig.7 Response surface for average wind speed of outlet 3 of multi-air-duct sprayer 
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4.2.2 出 风口 4 
对 出 风口 4 处 的 平均 模拟 风速 值 进行 分 析 ， 
得 到 回归 方程 如 公式 (4). 
Y,=33.02 +4.04 X 4-4.37 X B-1.31 X C- 
2.24 X 4B+3.78 X AC+1.72 X BC+ 
3.10 X Æ -2.21 X B’-2.28 X C? (4) 


将 回归 方程 的 显著 性 及 方差 分 析 置 于 
表 4 中 。 

由 表 4 可 知 ， 模 型 已 值 <0.01， 为 极 显著 ， 决 
定 系 数 R=0.9704， 表 明 模 型 拟 合 效果 较 好 ， 可 
用 来 分 析 及 预测 响应 值 。 该 模型 的 一 次 项 4、B 
REE, CH; 3 个 因素 之 间 两 两 交互 4C 极 显 
著 ，4B、BC 显 著 ; 二 次 项 4 极 显 著 ，B*?、C' 显 
著 。 出 风口 4 处 的 平均 风速 值 六 作为 评价 标准 ， 
在 数值 分 析 中 作为 优化 目标 ， 期 望 卫 最 大 化 ， 约 
束 是 导 流 板 1、 导 流 板 2、 导 流 板 3 的 长 度 ， 数 学 
模型 表示 如 公式 (5)。 

max: Y,=f(4,B,C) 


175 < A < 225 (5) 
St.) 35 <B< 85 
35<C< 60 


图 8 为 出 风口 4 平均 风速 值 的 响应 面 ， 可 以 
直观 地 发 现 ， 导 流 板 1、 导 流 板 2、 导 流 板 3 之 
间 的 交互 作用 明显 。 分 析出 风口 4 处 的 平均 风速 
值 ， 在 各 试验 因素 的 取 值 范围 内 ， 求 得 最 优 的 导 
流 板 长 度 寸 组 合 为 导 流 板 1 长 度 223.72 mm 、 导 
流 板 2 长 度 41.09 mm， 导 流 板 3 长 度 50.21 mm, 


ff Gms!) 


模拟 风速 


i 


(a) 导 流 板 1 与 导 流 板 2 交互 响应 面 


(b) 导 流 板 1 与 导 流 板 3 交 互 响应 面 


表 4 多 风 道 喷雾 机 出 风口 4 平均 风速 值 响应 面 回 归 方程 
的 显著 性 及 方差 分 析 结 果 
Table 4 Significance and variance analysis of response sur- 


face regression equation for outlet 4 of multi-air-duct sprayer 


方差 来 源 ”平方 和 IF FE P 值 ”显著 性 


模型 465.72 
A- 导 流 板 1 长 度 130.42 1 130.42 64.36 <0.0001 极 显 著 
B- 导 流 板 2 长 度 153.08 1 153.08 75.54 <0.0001 极 显 著 
C- 导 流 板 3 长 度 13.83 1 13.83 6.82 0.0348 显著 


51.75 25.54 0.0002 显著 


O 


AB 20.09 1 20.09 9.91 0.0162 显著 
AC 57.22 1 57.22 28.24 0.0011 极 显著 
BC 11.86 1 11.86 5.85 0.0461 显著 
A 40.49 1 40.49 19.98 0.0029 极 显著 
B? 20.51 1 20.51 10.12 0.0155 显著 
C 21.92 1 21.92 10.82 0.0133 显著 
残 差 14.18 7 2.03 
失 拟 项 14.18 3 4.73 
纯 误差 0.00 4 0.00 


总 离 差 479.90 16 


出 风口 4 平均 风速 值 为 43.310 m/s。 

为 检验 试验 结果 的 可 靠 性 ， 采 用 最 优 导 流 
板 参数 组 合 进行 3 次 平行 仿真 试验 ， 得 到 的 出 
风口 4 平均 风速 值 分 别 为 43.215、43.269 和 
43.478 m/s， 平 均值 为 43.321 m/s， 结 果 与 预测 值 
接近 ， 相 对 误差 为 0.02%， 表 明 该 模型 能 够 准确 
的 预测 仿真 结果 。 


(c) 导 流 板 2 与 导 流 板 3 交互 响应 面 
图 8 多 风 道 喷雾 机 出 风口 4 平均 风速 值 的 响应 面 


Fig. 8 Response surface for average wind speed of outlet 4 of multi-air-duct sprayer 


4.2.3 出 风口 5 
对 出 风口 5 处 的 平均 模拟 风速 值 进行 分 析 ， 


得 到 回归 方程 如 公式 (6)。 
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Y,=34.33-4.71 X A+3.54 X B-0.87 X C+ 
3.82 X AB+0.88 X AC+0.87 X BC- 
0.70 X A?- 1.69 X B?+0.69 X C? (6) 


将 回归 方程 的 显著 性 及 方差 分 析 结 果 列 于 
表 5 中 。 

由 表 4 可 知 ， 模 型 P 值 <0.01， 为 极 显 著 ， 决 
定 系 数 RR=0.9946， 表 明 模 型 拟 合 效果 较 好 ， 可 


RS 多 风 道 喷雾 机 出 风口 5 平均 风速 值 响应 面 回归 方程 
的 显著 性 及 方差 分 析 结 果 
Table 5 Significance and variance analysis of response sur- 


face regression equation for outlet 5 of multi-air-duct sprayer 


方差 来 源 ”平方 和 ABE HF FE PE 显著 性 


用 来 分 析 及 预测 响应 值 。 该 模型 的 一 次 项 4、B、 
C 极 显著 ; 3 个 因素 之 间 两 两 交互 48B 极 显著 ， 
4C、BC 显 著 ; 二 次 项 BF 极 显 著 ，4*?、C 显 著 。 
出 风口 5 处 的 平均 风速 值 关 作为 评价 标准 ， 在 数 
值 分 析 中 作为 优化 目标 ， 期 望 交 最 大 化 ， 约 束 是 
导 流 板 1、 导 流 板 2、 导 流 板 3 的 长 度 ， 数 学 模型 
表示 如 公式 (7)。 
max: Y,=f(4,B,C) 


175 < A < 225 (7) 
St.y 35 <B< 85 
35 <C< 60 


图 9 为 出 风口 5 平均 风速 值 的 响应 面 ， 可 以 


模型 363.88 9 40.43 143.45 < 0.0001 显著 
4- 导 流 板 1 长 度 177.14 1 177.14 628.50 < 0.0001 极 显著 
B- 导 流 板 2 长 度 100.46 1 100.46 356.43 < 0.0001 极 显 著 


KF 


KE 


C- 导 流 板 3 长 度 6.00 1 6.00 21.28 0.0024 极 显著 
AB 58.26 1 ”58.26 206.70 <0.0001 极 显著 
AC 3.09 1 3.09 10.96 0.0129 显著 
BC 3.03 1 3.03 10.75 0.0135 ”显著 
a 2.04 1 2.04 7.25 0.0310 ”显著 
B 11.99 1 11.99 42.55 0.0003 极 显著 
C 1.99 1 1.99 7.06 0.0326 显著 
残 差 1.97 7 0.29 
失 拟 项 1.97 3 0.66 
纯 误 差 0.00 4 0.00 


总 离 差 365.85 16 


HEL RUE Ges!) 


D 
a C: 导 流 板 3/mm 
B: 导 流 极 2/mm 


A: Fimm 


(a) 导 流 板 1 与 导 流 板 2 交 互 响应 面 


b y5 


(b) 导 流 板 1 与 导 流 板 3 交互 响应 面 


直观 地 发 现 ， 导 流 板 1、 导 流 板 2、 导 流 板 3 之 间 
存在 交互 作用 ， 且 交互 作用 明显 。 分 析出 风口 $ 
处 的 平均 风速 值 ， 在 各 试验 因素 的 取 值 范围 内 ， 
求 得 最 优 的 导 流 板 长 度 长 度 组 合 为 导 流 板 1 
长 度 175.00 mm 、 导 流 板 2 长 度 51.53 mm, & 
流 板 3 长 度 35.00 mm， 出 风口 5 平均 风速 值 为 
40.963 m/s. 

为 检验 试验 结果 的 可 靠 性 ， 采 用 最 优 导 流 
板 参数 组 合 进行 3 次 平行 仿真 试验 ， 得 到 出 风 
口 5 平 均 风 速 值 为 分 别 为 40.129、41.246 和 
41.875 m/s， 平 均值 为 41.083 m/s， 结 果 与 预测 值 
接近 ， 相 对 误差 为 0.29%， 表 明 该 模型 能 够 准确 
的 预测 仿真 结果 。 


PURE (ms!) 


D 
C: 导 流 板 3/mm 


A: 导 流 板 1/mm 


(c) 导 流 板 2 与 导 流 板 3 交 互 响应 面 


图 9 多 风 道 喷雾 机 出 风口 5 平均 风速 值 的 响应 面 


Fig. 9 Response surface for average wind speed of outlet 5 of multi-air-duct sprayer 


4.2.4 综合 分 析 
综 上 所 述 ， 由 于 出 风口 4 与 出 风口 3， 出 风 
口 3 与 出 风口 6 为 左右 对 称 出 口 ， 需 满足 对 称 位 


置 出 风 均 匀 。 为 减少 误差 、 确 保 以 上 数据 处 理 符 
合 要 求 。 确 保 试 验 结 果 的 准确 性 、 可 靠 性 ， 进 行 
相对 误差 分 析 ， 通 过 计算 以 上 三 组 参数 组 合 对 称 
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位 置 出 风口 风速 的 相对 误差 得 出 表 6。 由 表 6 可 FE 60.00 mm、 导 流 板 3 长 度 50.00 mm 为 最 优 参 
知 ， 确定 导 流 板 1 长 度 200.00 mm、 导 流 板 2 长 数组 合 。 
表 6 多 风 道 喷雾 机 对 称 出 风口 风速 相对 误差 分 析 


Table 6 Relative error analysis of symmetrical outlet wind speed of multi-air-duct sprayer 


组 别 出 风口 3/m:s” ”出 风口 6/m:s” ”相对 误差 /% MAKO sms 出 风口 5/m:s” ”相对 误差 /% 
200+60+50 39.135 41.320 5.58 33.022 34.328 3.95 
223.72+41.09+50.21 35.692 41.033 14.96 43.310 22.680 47.63 
175+51.53+35 38.141 40.301 5.66 35.953 40.963 13.94 


5 室内 风速 验证 试验 


5.1 试验 方案 


为 更 好 地 考察 果园 多 风 道 喷雾 机 送 风 系统 凡 \ 
场 的 分 布 情况 ， 以 仿真 优化 分 析 获 得 的 喷雾 参数 
组 合作 为 试验 条 件 ， 通 过 室内 风速 试验 对 送 风 系 
统 的 风 场 分 布 情况 进行 研究 ， 同 时 验证 仿真 模拟 
优化 结果 的 准确 性 。 

试验 仪器 包括 多 风 道 喷雾 机 送 风 系统 装置 、 
风速 测量 仪 (AS856$S， 和 硕 玛 公司 ) 、 标 杆 、 盒 
尺 、 三 项 异步 电动 机 (Y2-90L-2 极 ) 和 变频 器 
(M2-2R2G3) 等 。 将 标杆 距 地 面 0.50 m 处 作为 
第 一 个 监测 点 ， 以 此 为 基点 在 标杆 上 方 每 隔 
0.30 mm 进行 标记 作为 测量 点 ， 最 高 测量 到 3.50 m 
处 。 连 接 变 频 吉 与 三 项 异步 电动 机 ， 并 把 变频 器 
与 三 项 电源 相 接 ， 将 标杆 置 于 喷雾 机 出 风口 轴线 
1.25 mm 处， 打开 电源 ， 调 节 变 频 器 频率 由 0 逐渐 
增加 到 39 Hz， 此 时 风机 转速 为 2184 rmin， 待 风 
机 转速 稳定 后 打开 风速 测量 仪 ， 置 于 测量 点 处 对 
风速 进行 测量 ， 并 进行 记录 。 一 侧 测 量 完 成 后 对 
另 一 侧 按 同 样 的 测量 方法 进行 风速 检测 ， 每 侧 风 
速 分 别 做 三 次 平行 试验 。 


5.2 验证 结果 分 析 


记录 整理 喷雾 机 左右 两 侧 采 样 得 到 的 风速 
值 ， 单 侧 平 行 试验 后 的 风速 值 取 其 平均 数 ， 应 用 
数据 分 析 软 件 进行 处 理 得 到 如 图 10 的 点 线 图 。 
点 线 图 中 展现 了 在 距离 喷雾 机 出 风口 1.25 m 处 风 
场 的 分 布 情况 ， 以 高 度 为 纵 坐 标 可 以 发 现在 没有 
受到 果树 冠 层 干扰 的 情况 下 ， 风 速 由 上 层 到 下 层 


逐渐 增 大 ， 实 现 风 场 按 果 树冠 层 形状 分 布 。 依 据 
戴 奋 奋 六 所 述 未 速度 原则 ， 到 达 果 树 内 膛 的 气 
流风 速 应 达到 9 一 10 m/s， 在 室内 风速 实验 中 
上 层 平 均 风 速 为 15.750 m/s， 中 层 平均 风速 
为 20.830 mV/s， 下 层 平均 风速 为 28.270 m/s, W 
足 末 速度 原则 ， 并 且 经 试验 验证 喷雾 机 左右 两 侧 
风 场 分 布 均匀 对 称 。 


4.0 


一 一 左 侧 
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图 10 室内 风速 试验 风速 分 布 示意 图 
Fig. 10 Wind speed distribution diagram of indoor wind 
Speed test 


6 & it 


本 研究 针对 多 风 道 喷雾 机 内 部 气流 分 布 不 
均 ， 导 致 由 出 风口 吹出 的 气流 亲 乱 ， 不 能 使 雾 滴 
在 果树 冠 层 上 均匀 沉积 的 问题 SE Star-CCM+ 
软件 对 多 风 道 喷雾 机 内 部 气流 场 进行 仿真 分 析 ， 
确定 导 流 板 最 优 长 度 、 组 合 ， 并 进行 了 室内 风速 
试验 验证 。 

(1) 应 用 CFD 技 术 ， 基 于 StarCCM+ 软 件 对 
果园 多 风 道 喷雾 机 送 风 系 统 内 部 气流 进行 模拟 分 
析 ， 得 到 出 风口 1 一 6 的 风速 在 不 同 导 流 板 长 度 
的 标准 差分 别 为 0.7468、0.6776、 1.4441, 
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5.1305, 4.5768 Ail 0.8209. 

(2) 对 风速 标准 差 较 大 的 出 风口 3~S 进行 响 

面 分 析 ， 最 终 确 定 导 流 板 1 长 度 200.00 mm, 
导 流 板 2 长 度 60.00 mm、 导 流 板 3 长 度 $0.00 mm 
为 最 优 参数 组 合 ， 得 到 对 称 出 风口 3 和 出 风口 6 
的 风速 值 为 39.135 和 41.320 m/s， 相 对 偏差 为 
5.58%， 出 风口 4 和 出 风口 5 的 风速 值 为 33.022 和 
34.328 m/s， 相 对 偏差 为 3.95%， 最 终结 果 符 合 
设计 要 求 。 

(3) 通过 室内 风速 试验 表明 在 距离 喷雾 机 出 
风口 1.25 m 处 风 场 风速 由 上 层 到 下 层 逐 渐 增 大 ， 
实现 风 场 按 果 树冠 层 形状 分 布 ， 并 且 经 试验 验证 
喷雾 机 左右 两 侧 风 场 对 称 分 布 。 


参考 文献 : 

[1] 傅 锡 敏 , AE, 丁 为 民 . 我 国 果 园 植保 机 械 现 状 与 
技术 需求 [J]. 新 疆 农机 化 , 2011(1): 61-63. 
FU X, LYU X, DING W. Current situation and techni- 
cal demand of orchard plant protection machinery in 


China[J]. Xinjiang Agricultural Mechanization, 2011 
(1): 61-63. 

[2] BADULES J, VIDAL M, BONE A, et al. Comparative 
study of CFD models of the air flow produced by an 
air-assisted sprayer adapted to the crop geometry[J]. 
Computers and Electronics in Agriculture, 2018, 149: 
166-174. 

[3] 何 雄 奎 . 中 国 精准 施 药 技 术 和 装备 研究 现状 及 发 展 
建议 中. 智慧 农业 (中 英文 ) 2020, 2(1): 133-146. 


HE X. Research progress and developmental recom- 


mendations on precision spraying technology and 
equipment in China[J]. Smart Agriculture, 2020, 2(1): 
133-146. 

[4] 张 燕 妮 , HKI, 赵 娟 . 果园 风 场 与 雾 场 测量 建 模 方 
FELT]. 农机 化 研究 , 2020, 42(4): 264-268. 
ZHANG Y, ZHAI C, ZHAO J. Modeling and measur- 
ing method of orchard airflow and droplet field[J]. 


Journal of Agricultural Mechanization Research, 2020, 
42(4): 264-268. 

[5] PFEIFFER S A, GUEVARA J, CHEEIN F A, et al. Me- 
chatronic terrestrial LIDAR for canopy porosity and 
crown surface estimation[J]. Computers and Electron- 
ics in Agriculture, 2018, 146: 104-113. 

[6] GILE, ARNO J, LLORENS J, et al. Advanced technol- 
ogies for the improvement of spray application tech- 
niques in Spanish viticulture: An overview[J]. Sensors, 
2014, 14(1): 691-708. 

[7] PETROVIC D. Odnos selektivnog i konvencionalnog 
raspršivanja te njihov utjecaj na depozit i zanošenje 


[8] 


[9] 


[10] 


[11] 


[12] 


[13] 


[14] 


[15] 


[16] 


[17] 


[18] 


teku¢ine[D]. Osijek: Josip Juraj Strossmayer Universi- 
ty of Osijek, 2018. 

GRELLA M, MARUCCO P, MANZONE M, et al. Ef- 
fect of sprayer settings on spray drift during pesticide 
application in poplar plantations (Populus spp.)[J]. Sci- 
ence of the Total Environment, 2017, 578: 427-439. 
康 峰 , RIRA, 王 亚 雄 , 等 . 农药 雾 滴 沉 积 特性 研究 进 
展 与 展望 四. 农业 工程 学 报 , 2021, 37(20): 1-14. 
KANG F, WU X, WANG Y, et al. Research progress 
and prospect of pesticide droplet deposition characteris- 
tics[J]. Transactions of the CSAE, 2021, 37(20): 1-14. 
WU B, HUANG K, QING D, et al. Numerical simula- 
tion and optimization of perforated tube trolley in cir- 


cular cooler[J]. World Journal of Engineering and Tech- 

nology, 2017, 5(4): 684-695. 

ip, 丁 为 民 , Ehm, 等 . 3WZ-700 型 自 走 式 果园 

风 送 定向 喷雾 机 四 . 农业 机 械 学 报 , 2012, 43(4): 26- 

30, 44. 

QIU W, DING W, WANG X, et al. 3WZ-700 self-pro- 

pelled air-blowing orchard sprayer[J]. Transactions of 

the CSAM, 2012, 43(4): 26-30, 44. 

周 良 富 , 傅 锡 敏 , TAR, 等 .组 合 圆 盘 式 果 园 风 送 喷 

雾 机 设计 与 试验 四 .农业 工程 学 报 , 2015, 31(10): 

64-71. 

ZHOU L, FU X, DING W, et al. Design and experi- 

ment of combined disc air-assisted orchard sprayer[J]. 

Transactions of the CSAE, 2015, 31(10): 64-71. 

陈 帮 .离心 风机 进 气 箱 导 流 板 的 优化 设计 上 四. 煤矿 机 

械 , 2022, 43(2):130-134. 

丁 天 航 ， 曹 昌明 , 薛 新 宇 , SE. 果园 喷雾 机 单 双 风机 风 
道 气流 场 仿真 与 试验 四. 农业 工程 学 报 , 2016, 32 

(14): 62-68, 315. 

DING T, CAO S, XUE X, et al. Simulation and experi- 

ment on single-channel and double-channel airflow 


field of orchard sprayer[J]. Transactions of the CSAE, 
2016, 32(14): 62-68, 315. 

ENDALEW A M, DEBAERB C, RUTTEN N, et al. A 
new integrated CFD modelling approach towards air- 
assisted orchard spraying—Part I: Model development 
and effect of wind speed and direction on sprayer air- 
flow[J]. Computers and Electronics in Agriculture, 
2010, 71(1): 128-136. 

NUYTTENS D, ZWERTVAEGHER I K A, DEKEY- 
SER D. Spray drift assessment of different application 
techniques using a drift test bench and comparison 
with other assessment methods[J]. Biosystems Engi- 
neering, 2017, 154: 14-24. 

BADULES J, VIDAL M, BONE A, et al. Comparative 
study of CFD models of the air flow produced by an 
air-assisted sprayer adapted to the crop geometry[J]. 
Computers and Electronics in Agriculture, 2018, 149: 
166-174. 

HONG S W, ZHAO L, ZHU H. CFD simulation of air- 
flow in side tree canopies discharged from air-assisted 


ChinaXiv 合 作 期 刊 


85 


202302.00147v1 


chinaXiv 


sprayers[J]. Computers and Electronics in Agriculture, [22] 李建平 , 边 永 亮 , 杨欣 , 等 . 果园 多 风机 风 送 喷雾 机 作 


2017, 149: 121-132. 业 参 数 优化 与 试验 [J]. 吉林 大 学 学 报 (工学 版 ), 2022, 
[19] ZHENG Y YANG S, LIU X, et al. The computational 52(10): 2474-2485. 

fluid dynamic modeling of downwash flow field for a LI J, BIAN Y, YANG X, et al. Operational parameter 

six-rotor UAV[J]. Frontiers of Agricultural Science and optimization and testing of an air-assisted multi-fan or- 

Engineering, 2018, 5(2): 159-167. chard sprayer[J]. Journal of Jilin University (Engineer- 
[20] HOLOWNICKI R, DORUCHOWSKI G, ing Edition), 2022, 52(10): 2474-2485. 

SWIECHOWSKI W, et al. Variable air assistance sys- [23] 姜 宗 月 . 果园 定向 仿 形 弥 雾 机 的 研制 与 试验 [D]. 泰 

tem for orchard sprayers: Concept, design and prelimi- 安 : 山东 农业 大 学 , 2014. 

nary testing[J]. Biosystems Engineering, 2017, 163: JIANG Z. Development and experiment of directional 

134-149. profiling orchard mist sprayer[D]. Taian: Shandong Ag- 
[21] Kit, 张 燕 妮 , SEAS, 等 . 果园 风 送 喷雾 机 出 风口 ricultural University, 2014. 

K CFD 建 模 与 试验 加. 智慧 农业 (中 英文 )2021, 3 [24] RAA. 风 送 喷 雾 机 风量 的 选择 与 计算 四 . 植物 保 

(3): 70-81. 护 , 2008(6): 124-127. 

ZHAI C, ZHANG Y, DOU H, et al. CFD modeling and DAI F. Selection and calculation of the blowing rate of 

experiment of airflow at the air outlet of orchard air-as- air-assisted sprayers[J]. Plant Protection, 2008(6): 

sisted sprayer[J]. Smart Agriculture, 2021, 3(3): 70-81. 124-127. 


Design Optimization and Test of Air Supply System for 
Multi-Duct Sprayer 


GUO Jiangpeng, WANG Pengfei’, LI Xinhao, YANG Xin, LI Jianping, 
BIAN Yongliang, XUE Chunlin 


(College of Mechanical and Electrical Engineering, Agricultural University of Hebei, Baoding 071000, China) 


Abstract: In view of the uneven distribution of airflow inside the multi-air-duct sprayer, the air flow caused by the air outlet is 
disturbed and the droplet can not be evenly deposited on the fruit tree canopy. In this research, the length parameter of the inner 
baffle plate of the multi-duct sprayer was optimized. The Computational Fluid Dynamics (CFD) was used to simulate and ana- 
lyze the internal airflow of the air supply system of the multi-duct sprayer based on Star-CCM+ software. The standard devia- 
tions of the wind speed of the wind outlet 1~6 at different guide plates were 0.7468, 0.6776, 1.4441, 5.1305, 4.5768 and 0.8209, 
respectively. Among them, the standard deviations of wind speed value at Point 1, Point 2 and Point 6 were less than 1, indicat- 
ing that the change of deflector length has little impact on the speed change. The standard deviations of wind speed value at 
Point 3, Point 4 and Point 5 were large, indicating that with the change of deflector length, the wind speed at Air outlet 3, Air 
outlet 4, Air outlet 5 were greatly affected. On this basis, through the response surface analysis of Air outlet 3, Air outlet 4 and 
Air outlet 5, it was determined that, the length of Deflector 1 as 200 mm, the length of Deflector 2 as 60 mm and the length 
of Deflector 3 as 50 mm, was the optimal parameter combination. Under the optimal combination parameters, the wind speed 
values of symmetrical Air outlet 3 and Air outlet 6 were 39.135 and 41.320 m/s, respectively, with a relative deviations of 
5.58%. The wind speed values of air outlet 4 and air outlet 5 were 33.022 and 34.328 m/s, respectively, with a relative deviation 
of 3.95%, which meeting the design requirements of sprayer. The indoor wind speed test results showed that the average wind 
speed of the upper layer was 15.75 m/s, the average wind speed of the middle layer was 20.83 m/s, and the average wind speed 
of the lower layer was 28.27 m/s, which met the end speed principle. The wind field was distributed according to the shape of 
the fruit tree canopy. The wind field of the left and right sides of the sprayer was symmetrical distributed and the air distribution 
was uniform. The work can provide a reference for the design of multi-duct sprayer. 

Key words: Computational Fluid Dynamics (CFD); multi-duct sprayer; air supply system; flow field simulation; response sur- 


face methodology; uniform deposition 
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